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We report detection of giant cyclonic vortices in the gaseous disc of the spiral galaxy 
NGC 3631 in the reference frame rotating with the spiral pattern. A presence of such 
structures was predicted by the authors for galaxies, where the radial gradient of 
the perturbed velocity exceeds that of the rotational velocity. This situation really 
takes place in NGC 3631. 
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1 Introduction 



In previous papers [1] , [2] authors have shown that besides the spiral arms, well 
known for more than one and a half century in spiral galaxies as structures in 
the brightness distribution, there also exist structures, revealing themselves in 
the velocity field. The latter structures have the appearance of giant anticy- 
clonic vortices: their rotation is opposite to the rotation of the galactic disc. 
The centers of these vortices are located near the zero points of the force field 
(Lagrange points L 4 and L 5 , see [3]). In spiral galaxies these points are fixed in 
the close vicinity of the corotation circle, where the velocity of the rigid-body 
rotating spiral pattern coincides with the velocity of the differentially rotating 
disc. Being in such a position, the centers of the anticyclones turn out to be 
fixed with respect to the spiral arms and stationary if the spiral structure 
is stationary. If the spiral structure develops due to some instability in the 
disc, then the arms and vortices arise simultaneously and grow with the same 
growth rate [4]. As a result a unified spiral-vortex structure forms. 
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The growth of the perturbation amplitude leads rather often to the saturation 
of the instability, after which a spiral-vortex structure becomes stationary. 
From the very beginning this scenario implies that spiral arms in discs of 
galaxies have a wave nature, that satisfies modern theoretical conceptions [5] 
and agrees well with the results of the velocity field analysis in spiral galaxies 
([2], [6]). 

A typical rotation velocity field of a galactic disc is represented schematically in 
Fig. la. There is no way to close the trajectories of particles in the laboratory 
(inertial) reference frame by any radial component of the velocity (arrows 
radially directed in Fig. la). The situation is opposite in the rotating reference 
frame, where anticyclonic motion takes place under condition of appropriate 
azimuthal variation of the radial component (see Fig. lb). Note that in the 
latter case it is not necessary to have large radial velocities to form a vortex 
near the circle corotating with the reference frame due to small values of the 
azimuthal velocities in this region. 

It is not accidental that the radial velocity was chosen as a periodic function 
of the azimuth: in the spiral density wave such periodicity should take place 
both for the radial V r and for the azimuthal perturbed velocities. Arrows 
for are not shown in Fig.l. Taking into account the perturbed azimuthal 
velocity will not change the picture qualitatively, if 

\dV^/dr\ < \dV circ /dr\, (1) 

where Vd rc = V ro t — £l r f r is the circular velocity in the reference frame rotating 
with the angular velocity fl r f. Note that even a weak inequality is quite enough 
to form only anticyclonic vortices like in Fig. lb (for more details see [7]). 
Under the opposite condition: 

\dV+/dr\ > \dV circ /dr\ (2) 

not only anticyclonic, but also cyclonic vortices can appear [8]. In the next 
section we will consider both these cases in a schematic model of a galactic 
disc with a well-defined two-armed spiral structure (Grand-Design galaxy). 



2 Qualitative description of anticyclone and cyclone formation in 
a disc of Grand-Design galaxy 

Let us assume that two-armed Grand-Design spiral structure exists in a dif- 
ferentially rotating gaseous disc. Then the functions describing the perturbed 
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surface density and velocity components, related to the density wave, may be 
approximated in the form: 

a(r, p) = C a (r) cos[2y9 - F a (r)} , (3) 

V r (r, <p) = C r (r) cos[2v? - F r (r)} , (4) 

V v (r,<p)=Cy(r)cos[2<p-Fy(r)]. (5) 

As we consider the momentary picture of the spiral structure, a time depen- 
dence is absent in the equations above. 

It is natural to ask, whether it is possible to use expressions (3)-(5) for the 
description of the dynamics of spiral density waves with a finite amplitude. 
We have checked it in [6], [14] by the example of the galaxy NGC 3631, where 
cyclones have been found just now (see Section III below). The operating 
sequence of the check was the following. First, we analysed the azimuthal ex- 
pansion in the harmonic series of the observed surface brightness of this galaxy 
at different galactocentric radii. A domination of the second Fourier harmonic 
over others (see Fig. 1 in [14]) proved the correctness of the approximation 
expressed in Eq.(3). 

Then the observed line-of-sight velocity field V obs (r, ip) was analysed. This 
velocity can be expressed as (see [2], [9]): 

V obs (r,p) — V s + Vtp(r,ip) cosy? + V r (r,ip) simp sini + V z {r, ip) cos«.(6) 

Here i is an inclination angle (angle between the rotation axis of galaxy and 
line of sight) and V s is systemic velocity of a galaxy (velocity of the center of 
mass of the galaxy). 

Within the frame of the present model 

V r (r,<p) = V r (r,(p), 
V v (r,<p) = V rot (r) + V v (r,<p), (7) 
V z (r,<p) = V z (r,p) = C z (r) cos[2p - F z {r)\ , 

where V r , V v are determined by (4) and (5). Substituting (7) in (6) one obtains 
the model representation of the Fourier expansion for V obs (r, p). 

Equating the coefficients of the model expansion to the coefficients calculated 
from observational data, 

V obs (r,p) = V s + £(af cosr^ + fcf sinrup) , (8) 
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we obtain the following set of relations between unknown characteristics of the 
vector velocity field of the galaxy and the Fourier coefficients of the observed 
line-of-sight velocity field: 



V rot + - (C r sin F r + C v cos F v ) 




(9) 



— C r cos F r + Cr, sin F, 



2b' 



>l bs , 



(10) 



C z cosF z cosi = a° 2 , 



(11) 



C z sin F z cos i 




(12) 



— C r sin F r + Cu, cos F, 



2a' 



3 > 



(13) 



C r cosF r + C^sinF^ = 2bf s . 



(14) 



From (9)-(14) we see that the perturbed velocity components V r ,Vy and 14, 
written in the form of (4), (5) and (7) give contributions only in the first three 
harmonics of the Fourier expansion. This is really observed in the case of NGC 
3631 (see Fig. 3 in [6]), where the first three Fourier harmonics of the line- 
of-sight velocity field dominate over the higher-order harmonics. This gives 
evidence for the correctness of the representation of the perturbed velocity 
components in the form (4), (5) and (7) for the galaxy NGC 3631. 

To draw schematically the spiral-vortex structure in different reference frames 
we need relations between the phases F r (r), F v (r) and F a (r). These relations 
were derived in [9] from the set of linearized equations. In spite of this sim- 
plification, an excellent qualitative agreement was found between the phase 
of the "modified" third harmonic (see below) of the observed line-of-sight 
velocity field and F a (r) (see Fig. 5 in [6]), which enables to use the phase 
relations obtained in [9] for a qualitative description of the residual velocities. 
The modified third harmonic is a two-armed spiral with phase F 3 — tt/2, i.e. 
oc cos(2(/9 — F 3 + 7r/2), where F 3 is the phase of the original radial-velocity-field 
third harmonic. 

Below a list of these relations is presented. 

a) Relations between the phases of the radial velocity V r and the surface 
density perturbation a. 

• On the corotation circle r = r c : 




(15) 



4 



whence 

V r (r c ,<p) = -C r (r c ) sin[2y9 - F a (r c )\, (16) 

i.e. on the corotation circle the extrema of V r coincide with the zeroes of a 
and vice versa. 

• Inside the corotation circle, r < r c , 

F r = F a + vr, (17) 

whence 

V r (r,<p) = -C r (r)cos[2ip-F a (r)\, r < r c , (18) 

i.e. inside the corotation circle the maxima of V r coincide with the minima 
of er and vice versa. 

• Outside the corotation circle, r > r c , 

F r = F a , (19) 

whence 

V r (r,<p) = C r (r)cos[2<p-F a (r)}, r > r c , (20) 

i.e. outside the corotation circle the maxima and minima of V r coincide, 
respectively, with the maxima and minima of a. 

b) Relation between the phases of the azimuthal velocity V v and the surface 
density a. 

• In all cases, both on the corotation circle, r = r c , and out of it, r ^ r c , we 
have 

F v = F v + |, (21) 

therefore 

V v (r,<p) = C^r)sm[2 V -F a (r)], (22) 

i.e. the extrema of coincide with the zeroes of a and vice versa throughout 
the disc. 

Using the expressions (16), (18), (20), and (22) we can draw a scheme for 
the velocity field with vortices overlayed on the scheme of the surface density 
distribution (Fig. 2). In Fig. 2b we can see the presence of cyclones and an- 
ticyclones in the field of residual velocity (that corresponds a reference frame 
locally corotating with the disc at each radius, V circ (r) = 0). This fact by 
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itself does not depend on the amplitude of the spiral density wave, although 
the sizes of the vortices depend on this amplitude. 

Lindblad and Langebartel [10] were the first to calculate the field of displace- 
ments of the stars in the gravitational potential of a bar. The form of the star 
displacement field resembles the system of two cyclones and two anticyclones, 
which is similar our schematic Fig. 2b. Refering on [11] Lynden- Bell [12] has 
noted the existence of cyclonic and anticyclonic trajectories in Fig. 9 of [10]. A 
similar picture of four vortices in the velocity field of a gaseous disc - two cy- 
clones on the bar and two anticyclones between the spiral arms - was obtained 
in a recent paper [13]. 

In the reference frame rotating with the angular velocity il ph of the two-armed 
spiral pattern the picture of vortices differs qualitatively from the previous one. 
In the case (1) when the gradient of the perturbed azimuthal velocity is less 
than the circular velocity gradient, i.e. the anticyclonic shear dominates over 
the cyclonic one, the appearance of cyclones is impossible - we can see only 
two anticyclones with centers near the corotation circle between two spiral 
arms (Fig. 2c). In the opposite case (2) when the gradient of the perturbed 
azimuthal velocity dominates over the circular velocity gradient, the location 
of the cyclones is determined by the place of this domination. Cyclones can 
survive almost in the same places as in Fig. 2b, i.e. on the corotation circle, 
but their sizes will be smaller. There are also two other variants when the 
cyclone center moves either inside or outside the corotation circle along the 
zero lines of the radial velocity. Finally, both last these variants can be realized 
simultaneously. In this case four cyclones exist in the velocity field as is shown 
in Fig. 2d. 



3 Discovery of Giant Cyclones in the Gaseous Disc of the Spiral 
Galaxy NGC 3631 with the 6m telescope in SAO 

Finding giant vortices in a spiral galaxy requires the reconstruction of the 
vector velocity field in the reference frame rotating with the angular velocity of 
the spiral pattern. Thus, first of all the corotation radius should be determined. 
In this paper for NGC 3631 we use the position of the corotation derived in 
[14] on the base of the Fourier analysis of the observed line-of-sight velocity 
field. The observations of NGC 3631 were carried out with the scanning Fabry- 
Perot interferometer with Russian 6m telescope. The details of observations 
were described in [6]. 

It is necessary to determine V r and V v = V ro t + V 9 , i.e. to find five unknown 
functions: V rot (r), C r (r), C v (r), F r (r), F^r) (see Eqs. (4), (5), and (7)). These 
five functions are connected with the characteristics of the observed velocity 
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field by the four relations (9), (10), (13), (14). An additional condition required 
to close the system should have a theoretical origin. Unfortunately, up to now a 
reliable condition valid for any density wave amplitude is not available. Several 
possibilities discussed in the literature [2], [9] have limited applicability. To 
overcome this difficulty, we propose the following approach. 

Among the functions listed above, V ro t(r) could be the most reliably estimated 
on the base of independent observational data. For this purpose the equilib- 
rium condition of the gaseous disc rotating in a gravitational potential ^ can 
be used 

= ^- (23) 
r or 



The latter is determined from the mass distribution in a galaxy or its sur- 
face brightness maps assuming a costant mass-to-light ratio to be known. To 
calculate the right hand side of the equation (23) we use a three-component 
dynamical model of a spiral galaxy similar to [15]. In spite of the roughness 
of the model, the V rot curve found by this method is in the region between 
(oi)min and (ai) ma x (see Fig. 3; a peculiarity of the thick curve is described 
below) . 

The same result is obtained by another way. From equation (9) it follows that 
the difference \a\ — V rot \ cannot exceed the amplitudes C r and C^, which, in 
turn, are connected by equations (10), (13) and (14) with the Fourier coeffi- 
cients 6° fcs , a^ bs , and b^ bs , determined from our observations. In Fig. 4 the radial 
behaviour of bf s — b obs = C r cos F r and bf s +bf s = sin F v is demonstrated. 
The extrema of these functions allow to estimate the amplitudes C r and C v . 
According to Fig. 4, one can conclude that in NGC 3631 a maximum value of 
the amplitude of the residual velocities can reach at 60 km/s, i.e. 

I^l Vrot | max 

< 60 km/s. (24) 



The conditions (23) and (24) do not allow to calculate the function V rot (r) 
exactly. Nevertheless, they set limits on the variations of both the amplitude 
and the form of V rot (r). Within these limits, we choose a set of trial curves 
(Fig. 3) and analyse the velocity field obtained from Eqs. (9)-(14) for a given 
V rot (r). 

The general views of the restored velocity fields are qualitatively similar for 
the full set of the rotation curves. In particular, in all cases, except the most 
extreme ones, distinct areas of anticyclonic as well as cyclonic shear flows exist. 
However, only rare examples demonstrate the presence of areas with trapped 
fluid particles moving along stream lines bounded by closed separatrices. Such 
behavior could be a natural consequence of the quasi-stationary character of 
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the density wave. It is an argument in favour of choosing just these examples of 
V ro t( r ) as the best-corresponding to the real dynamical model of the galaxy. 
The best example of U rot from the point of view of the stationarity of the 
stream lines is shown by the thick curve in Fig. 3. In other examples the vortices 
have not close separatrices and their stream lines are spirals - strongly inwards 
or outwards. 

In Fig. 5 we can see the amplitudes of radial (V r ) and azimuthal (V^) com- 
ponents of the residual velocity. These amplitudes, C r and were obtained 
from four equations (9), (10), (13), and (14) after substituting in them V rot (r) 
in the form drawn in Fig. 3 by the thickest line. 

Using the functions C r and C v in Fig. 5 let us show, that the necessary condi- 
tion (2) of the cyclone formation is satisfied in the galaxy NGC 3631 outside 
and inside the corotation circle and is not satisfied on the corotation circle. In 
the reference frame, rotating with the angular velocity of the spiral pattern 
(Q = Q p h), the inequality (2) can be rewritten in the form 





> 






d 


dr 


dr 




dr 



'rot 



tt ph r) 



(25) 



Substituting here the expression (5) for V v , we obtain: 

C'^coaptp-Fj + F^CySmptp-Fy) > \V^ C \ , (26) 



where "prime" denotes the derivative with respect to r. 

In order that cyclones should be possible the condition (26) should be fulfilled 
at least near the center of the vortex. In the centers 



V r 



V r = C r cos(2(/9 — F r 







(27) 



This equation gives a possibility to reduce the left hand side of inequality (26) 
to a form more suitable for an estimation. 

In the vicinity of the center of the cyclone in residual velocities (Fig. 2 b) on 
the corotation circle from (15) and (21) we have F v — F r + ir. Thus with the 
help of (27), the condition (26) is reduced to 



(28) 



As follows from Fig. 5 b at the corotation radius (r ps 42 arcsec) C v ~ 0. 
Hence, as the value of F' 9 is very small the inequality (28) cannot be fulfilled. 
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In the case of "external" cyclones (see Fig. 2d) from (19) and (21) we have 
F ip = F r + tt/2. Then using (27) from (26) we obtain 



> 



\V- I 

I arc I 



(29) 



This relation explains the coincidence of the location of the external cyclones 
in NGC 3631 (52 < r < 59 arcsec in Fig. 7) with the region of rapid growth 
of C,p (50 < r < 60 arcsec in Fig. 5b). 

In the case of the "internal" cyclones from (17) and (21) it follows F v = 
F r —7r/2. Then in this case the condition (26) also has the form (29). According 
to Fig. 5b it is fulfilled in the region 22 < r < 31 arcsec, which really contains 
the internal cyclone region 25 < r < 30 arcsec (Fig. 7). 

The vector velocity field of the galaxy NGC 3631, when V rot (r) is taken in the 
form presented in Fig. 3 by the thickest line, is shown in Figs. 6 and 7. In the 
figures some streamlines of the velocity fields are presented. If the velocity field 
is stationary the trajectories of fluid particles coincide with the correspondent 
streamlines. One can see that the residual velocity field of NGC 3631 (Fig. 6) 
corresponds to that of Fig. 2b and the full velocity field in the reference frame 
of the spirals (Fig. 7) corresponds to either of two cases presented in Fig. 2d 
when the centers of two cyclones lie outside the corotation circle. Internal 
cyclones in the case shown in Fig. 7 are not enveloped by close streamlines. 



4 Conclusions 



I. An analysis of velocity fields in Grand-Design galaxies shows that: 

1) the field of residual velocities contains two cyclones and two anticyclones 
with centers on the corotation circle; 

2) in the reference frame rotating with the spiral pattern the velocity field 
belongs to one of two types: 

a) under the condition (1) it contains only two anticyclones with centers near 
the corotation circle; 

b) under the condition (2) besides the anticyclones the field contains also 
either two or four cyclones. In the former case the cyclone centers lay either 
on, or outside, or inside the corotation circle. In the latter case, two pairs of 
cyclones appear with centers inside and outside the corotation circle. 

II. Our analysis of the velocity field data obtained by our team at the 6m 
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telescope in the Special Astrophysical Observatory of the Russian Academy 
of Sciences shows that the Grand-Design galaxy NGC 3631 belongs to the 
type b). 
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Figure Captions 



Figure 1. A scheme for the anticyclonic vortices formation in a reference frame 
rotating with an arbitrary angular velocity. Unperturbed rotation velocities 
and radial perturbed velocities are marked by solid arrows and dashed arrows 
correspondingly, a) In the laboratory system of coordinates the rotation ve- 
locity of the disc varies with radius without change of sign. In this case the 
radial perturbed velocity can not provide a vortex formation despite the reg- 
ular change of its sign along azimuth, b) In a rotating reference frame, where 
the dotted circle is at rest, the same field of radial perturbed velocities as in 
Fig. la participates in the anticyclone formation. 



Figure 2. A scheme for the velocity field with vortices in different reference 
frames. Vectors of unperturbed rotation velocity field are shown by solid ar- 
rows, radial and azimuthal components of the residual velocity field - - by 
dashed arrows. Solid curves of different thickness - the thickest, the thinnest, 
and intermediate - tracing the azimuthal locations of the maxima, minima 
and zero values of the perturbed surface density defined at every radius are 
denoted by a max , a m i n , and 0(5-) respectively. A and C denote anticyclones 
and cyclones correspondingly, a) The angular velocities of a disc (Q(r)) and 
of the spiral pattern (Q p h) in the laboratory reference frame. The dashed line 
represents the corotation circle, b) The residual velocity field, which is the 
result of the subtraction of the rotation velocity from the full velocity field. 
Dotted lines show the boundaries of vortices, c) The velocity field in the ref- 
erence frame rotating with the angular velocity Q p h of the two-armed spiral 
pattern. Only two anticyclones can be seen in the vicinity of the corotation 
circle, when the gradient of perturbed azimuthal velocity is lower than the 
rotation velocity gradient, d) The velocity field in the same reference frame 
as in the previous figure, but in the case, when the gradient of the perturbed 
azimuthal velocity exceeds the rotation velocity gradient. In this case one can 
see two cyclones and two anticyclones. Depending on the residual velocity field 
geometry, the domination of the cyclonic shear over the anticyclonic one can 
take place either near corotation, or inside (or outside) the corotation circle. 
Figure d) demonstrates one more possibility when four cyclones exist simul- 
taneously - two inside and two outside the corotation circle together with two 
anticyclones situated at corotation. 



Figure 3. Examples of trial curves used to represent the rotation curve (V ro t(r)) 
in NGC 3631 are shown by solid lines together with the observed behaviour of 
(triangles). The thickest line marks the rotation curve corresponding 
to a quasi-stationary gaseous disk. 
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Figure 4. The radial dependence of bf s — b° bs = C r cos F r and bf s + b° bs 
= sin F v observed in the spiral galaxy NGC 3631. An estimation of the 
amplitudes of the velocity components from the extrema of the presented 
functions gives max{C r ) ~ maxiC^) ~ 60 km/s. 

Figure 5. Some parameters of the velocity field which correspond to a quasi- 
stationary regime of the vortex structure. 

(a) The calculated amplitude of the radial velocity as a function of r. 

(b) The same for the azimuthal velocity (squares) with overlaid profile of 

\V cir c(r)\. 

Figure 6. The residual velocity field in the plane of the gaseous disc of the 
galaxy NGC 3631. One can see two cyclones and two anticyclones. Solid lines 
show some streamlines of the residual velocity field. 

Figure 7. The full velocity field of the same galaxy in the reference frame ro- 
tating with the angular velocity of the two-armed spiral pattern. Overlayed 
squares show the position of maxima of the second Fourier harmonics of bright- 
ness map of NGC 3631 in H Q line. Thin circle shows the position of the corota- 
tion. Presented streamlines (solid curves) were calculated using reverse "time" 
direction that allows to reveal separatrices evidently. Two cyclones and two 
anticyclones are enveloped by close streamlines. Other two cyclones are not 
enveloped by close streamlines. The anticyclone centers lie almost at corota- 
tion; the cyclone centers lie outside corotation Cyclones lie in the vicinity of 
spiral arms, anticyclones are situated between them. The agreement of the 
vortices position with theoretical predictions (scheme in Fig. 2d) is very good. 
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